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;-^, | Abstract 

^ ■ We show how a mass mixing matrix can be generated dynamically, for 

Q ■ two massless fermion flavours coupled to a Lorentz invariance violating (LIV) 

t^^ . gauge field. The LIV features play the role of a regulator for the gap equa- 

. [ tions, and the non-analytic dependence of the dynamical masses, as functions 

Q I of the gauge coupling, allows to consider the limit where the LIV gauge field 

cn ' eventually decouples from the fermions. Lorentz invariance is then recov- 

ered, to describe the oscillation between two free fermion flavours, and we 
check that the finite dynamical masses are the only effects of the original 
^ ! LIV theory. We also discuss briefly a connection of our results with the case 

j^ I of Majorana neutrinos in both, the standard model, where only left-handed 

(active) neutrinos are considered, and extensions thereof, with sterile right- 
handed neutrinos. 



1 Introduction and Motivation 

The generation of quark, lepton and vector boson masses, as described 
in the Standard Model due to their coupling with the Higgs boson (Spon- 
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taneous Symmetry Breaking), seems to have been confirmed by the latest 
experimental results at the Large Hadron Collider pQ, with the discovery 
of a Higgs-like (scalar) particle. However, the origin of neutrino masses is 
still not well established, although the seesaw mechanism seems the most 
elegant and simple for such a purpose [2] . Seesaw mechanisms involve neces- 
sarily Majorana mass type fermions, and heavy right-handed neutrino states, 
without standard model interactions (sterile), whose exchange explains the 
smallness of the active neutrino (left-handed) species of the Standard Model. 
Such sterile neutrinos have not yet been discovered in Nature [3]. 

The possibility, therefore, of generating neutrino masses dynamically with- 
out the involvement of heavy right-handed states is still at play. In this ar- 
ticle we take some preliminary steps in this direction and envisage scenarios 
in which flavour oscillations can arise dynamically, from the flavour-mixing 
interaction of two massless bare fermions with an Abelian gauge field, which 
has a Lorentz-Invariance- Violating (LIV) propagator. Lorentz symmetry vi- 
olation is achieved by higher order space derivatives, which are suppressed 
by a large mass scale M. This mass scale allows the dynamical generation 
of fermion masses, as was shown in |1] with the Schwinger- Dyson approach. 
Another role of this mass scale is to lead to a finite gap equation, and there- 
fore to regulate the model. Further studies using a similar model were done 
in [5] to generate a fermion mass hierarchy. 

Moreover, LIV U(l) gauge models of the form suggested in [1] have been 
shown [B] to arise in the low-energy limit of some consistent quantum gravity 
theories, for instance when the U(l) gauge theory is embedded in a stringy 
space time foam model, with the foamy structures being provided by (point- 
like) D-brane space-time defects ("D-particles"). In such microscopic models, 
the gauge field was one of the physical excitations on brane world universes 
interacting with the D-particles. It was observed in [B] that the LIV La- 
grangian of |1] can be obtained from a Born-Infeld-type Lagrangian of the 
U(l) gauge field in the D-particle background, upon an expansion in deriva- 
tives. Lorentz Violation arises locally in such models as a result of the recoil of 
the D-particle defects during their interaction with open strings representing 
the U(l) excitations. Other works involving quantization of higher-order- 
derivative extensions of Quantum Electrodynamics can be found in [7] . 

An important point is the following structure of the dynamical fermion 
mass |1] 

rudyn - M exp(-a/e^) , (1) 



where a is a positive constant and e is the gauge couphng. Such a non- 
analytical form is well-known in the studies of magnetic catalysis [S], and 
it can be derived from a non-perturbative approach only, as the Schwinger- 
Dyson derivation of a gap-equation, used in |3] and here. From the expression 
([1]), one can see that it is possible to take the simultaneous limits 

M -)- cx) and e -)> , (2) 

in such a way that the dynamical mass ([1]) remains finite, corresponding to a 
physical fermion mass. This procedure is consistent in the string-embedding 
case of [B], where Lorentz symmetry is recovered in the limit of vanishing 
density of D-particles. In that model, the LIV scale can diverge in the case 
of vanishing D-particle density and zero fluctuations of the recoil velocity 
(evaluated over a stochastic population of D-particles), in which limit also 
the coupling can go to zero [H], in such a way that the dynamically generated 
fermion mass remains finite. This is a physical case in which the vector U{1) 
fields appear as LIV regulators, implying dynamical mass for fermions. 

In the present article we shall consider the regularisation ([2]) in a more 
generic sense. We note at this point that the role of Lorentz symmetry 
breaking as an UV regulator of a quantum field theory has been considered 
in [inj, but from a rather different perspective than ours. Our aim here is 
to discuss the dynamically generated mass for fermions and/or the induced 
oscillations among fermion species, using the coupling of the fermions with 
such a LIV regulator gauge field. 

Oscillations of massless neutrinos were already studied in [11], where neu- 
trinos are considered open systems, interacting with an environment. Such 
oscillations have also been studied in [12], in the framework of LIV models, 
involving non-vanishing vacuum expectation values for vectors and tensors. 
Other constraints and consequences of these LIV models are given in [T3] . 
Whilst these studies have been questioned by phenomenological constraints 
[2], our present model, based on higher order space derivatives, is not ex- 
cluded. We note that dynamical generation of flavour oscillations was also 
studied in the context of Lifshitz theories [15] , and a detailed analysis of this 
mechanism was done in [16], for two Lifshitz fermions coupled by a renor- 
malizable four-fermion interaction. 

In the limit (^, the non-physical gauge field decouples from the theory, 
and hence the gauge dependence of the dynamical mass is avoided (although 
this problem can be understood perturbatively in the framework of the pinch 



technique [T7], as explained in [B]). We stress here an essential feature of 
the mechanism described in the present article. Although LIV operators 
are suppressed by a large mass scale, so that the corresponding effect is 
negligible at the classical level, quantum corrections completely change this 
picture, and lead to finite effects. In our present study, the finite effect is the 
dynamical generation of fermion masses, which is present even after setting 
the LIV-suppressing mass scale M to infinity. Note that the order of the steps 
followed is important: quantization is done for finite mass M and coupling 
e, after which the simultaneous limits ([2]) are taken. 

The structure of the article is the following: Next section [2] introduces 
the model and derives the corresponding gap equations which must be sat- 
isfied by the dynamical masses. We consider the corresponding constraints 
and calculate the dynamical masses in the relevant cases in section [31 In 
subsection 13.61 we discuss the "Lorentz- Invariant limit" (|2]), in which the LIV 
gauge field decouples from fermions, and we demonstrate that relativistic 
dispersion relations for fermions are indeed recovered. The extension of the 
Dirac fermion case to chiral Majorana fermions, as appropriate for neutrinos 
either in the standard model or in seesaw-type extensions thereof, involving 
sterile neutrinos, is discussed in section HI Finally conclusions and outlook 
are presented in section |5l Technical aspects of our work are given in two 
Appendices. 

2 Dynamical Fermion Mass Matrix 

2.1 The Field Theory Model 

The LIV model we consider is 

^ = -\fA^ - -^)F^'' + ^(^^ - r4)^, (3) 

where F^j, is the Abelian field strength for the gauge field A^ and A = —did^ 
is the Laplacian (the metric used throughout this work is diag(l, -1, -1, -1)). 
The mass scale M suppresses the LIV derivative operator A, and can be 
thought of as the Plank mass, which eventually will be set to infinity. \Ef is a 
massless fermion doublet 



and the flavour mixing matrix r features tlie gauge couplings (ei, 62, e) as 
fei -ie\ 61 + 62 ei -62 , . 

where ctj are the usual Pauli matrices and 1 is the 2x2 identity matrix. 
The fermions ipi and ip2 in (E]) are Dirac, but the structure of the gap equa- 
tions that will be derived bellow remains the same in the case of Majorana 
fermions, hence the corresponding dynamical masses are independent of the 
nature of fermions. As already noted in the previous section, the Lagrangian 
([3]) can be derived from a stringy space time foam model, as shown in [6]. 
We mention in passing that such a space time foam model was already used 
to study decoherence in flavour oscillations, both in flat space time and in a 
Friedman- Robertson- Walker metric |18j . 
The gauge field bare propagator is 

D = ^- 1^ ^^^ +C ^^^ ^ (6) 

where C is a gauge fixing parameter, which appears in the final expression 
for the dynamical masses, but do not play a role in the simultaneous limits 

M — )■ 00 and 61, 62, e — > , (7) 

that leave the dynamical masses finite, as we discuss further on. 

We note that the flavour mixing interaction \l/r^\Ef can be at the origin of 
a gauge boson mass, which is dynamically generated, as fermion masses. This 
alternative to the Higgs mechanism is explained in [12], and was extended 
to a LIV model in [20] • In the present article, we disregard the possibility 
to generate a gauge boson mass dynamically, since, as we shall demonstrate 
below, the flavour mixing coupling e vanishes necessarily for consistency of 
the model in the case there is dynamical generation of fermion oscillations. 

The bare fermion propagator is S* = ip/p^, where p^ = {00, p), and we 
assume the dynamical generation of the fermion mass matrix 

fmi n\ mi+m2 mi-m2 , , 



with eigenvalues 



rui 



+ m2 iy/(mi - 7722)^ + 4/i 



2 



Xm± = \ ± ^^ ,' . (9) 



Because the mass matrix contains in general non-diagonal elements, the 
flavour eigenstates I-j/^j), i = 1,2 are not the same as the mass eigenstates \ip±) 
and there is mixing and oscillations, provided the energy eigenvalues E± = 
yjjp- + A^ are different [J. As usual, the flavour eigenstates are connected to 
the mass (energy) eigenstates by a unitary transformation, parametrised by 
a mixing angle Q: 

' ^\\ _ ( cos9 sin6' ^ f'^+\ f-.r.\ 



^^. 



-sin6' cos^ 



^- 



and if at time t = one has a flavour 'il)i{t = 0) then the probability of 
obtaining (under Hamiltonian evolution) the other flavour ip2{t) at t > is 
non trivial and given by 



Vi2{t) =sin^2^sin^ 



E.-E^ 



t 



(11) 



and the survival probability Vu = 1 — Vu- These constitute the flavour 
oscillations. We stress that, as becomes evident from flTTl) . non-trivial mixing, 
6' 7^ 0, is not sufficient for oscillatory behaviour among flavours, one needs 
necessarily different energy levels E_^ ^ E_ as well. In what follows we shall 
identify cases where mixing and/or oscillations are generated dynamically, 
as a result of the coupling of the fermions with the LIV gauge bosons. 

If we neglect other quantum corrections, the dressed fermion propagator 
G, obtained by solving the equation 



-M) = il 



is then 



G 



P^ + f{mi + 7722) + ?7li?7l2 — /U^ 

- 2/i^(p^ + mim2) + /i^ 

mi — 7712 



(jp, _ ffi'^'^ijP' — 77^2) 



(12) 



(13) 



7711 + '^2n 



-0-3 + ^Cfi 



We must check in what follows that the dynamical masses ?7ii, 7712, /i assumed 
here can indeed be generated by quantum corrections, which is obtained using 
the Schwinger-Dyson approach. 



^We shall check in subsection 13.61 that the relativistic dispersion relations for the 
fermions are indeed obtained in the Lorentz Invariant Limit ([7]), after (finite) dynami- 
cal mass generation. 
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2.2 Schwinger-Dyson Gap equations 

The self-consistent Schwinger-Dyson equation for the fermion propagator 
has the usual structure [21], and is not modified by the LIV term in the 
Lagrangian ([3]). If we neglect corrections to the wave functions, the vertices 
and the gauge propagator, the Schwinger-Dyson equation reads for our model 



G 



-1 



S~' 



D^, tY G tY ■ 



(14) 



The previous loop integral is finite as a consequence of the LIV term p^/M"^ 
in the denominator of the gauge propagator ([H]). We show in Appendix A 
that the equation flTll) leads to the following four gap equations, which must 
be satisfied by the three masses mi,m2,/i, 



mi 

4 + C 

7722 

4 + C 

4 + C 




where 



(e^mi + e^m2)h + (/i^ - mim2){e\m2 + e^mi)l2 

{e\m2 + e^mi)Ii + (/i^ - ■mim2){elmi + e^m2)l2 

/i(eie2 - e^)[/i - {^^ - 71111712)12] 

e(eimi + e2m2)/i + e(/i^ - mim2)(eim2 + e2mi)/2 

J{Al) - J{Al) 



(15) 



h 
h 



Al 



^2 



(16) 



Al-A^_ 



J{At 

Al 



A^_ 



and 

J{Al] 

Al 



dp- 



ff 



47r3 Jq "^ 1 + j?/M^ 
2 I ^2 



duj 



u^ 



jp. 



co'^ + fp 



Al 



nil + 7712 + 2/^^ v^(m^ — rri^Y + Aii^{mi + 7722)^ 



(17) 



2 2 

After the integration over the frequency u and the momentum p in the inte 
grals J(A^), we obtain for M » 7711,7712, fi 

1 1 



h 



167r2 A\ - Al 
1 1 






167r2 Ai 



J\2 



(18) 



The four equations fllSp do not have obvious solutions, since they must be 
satisfied by only three unknowns mi, 7722, /i. In what follows, we study dif- 
ferent solutions. The ones allowing for the generation of fiavour oscillations 
must have /i 7^ 0. 

2.3 Constraints 

From the first two equations flT^ . one obtains for efe^ 7^ e^: 

r ^ 1 elml - elmj 

A + C {elel-e^){ml-ml) ^'^^ 



(/i^ - mim2)/2 



1 mim2\e'{ — 62) + e (mg — ^^i 



4 + C {elel-e^){ml-ml) 

and the third and forth equations lead to the following constraints respec- 
tively 

/i(mi +m2)(e2mi + eim2)(ei - 62) = (20) 

e(e2mi + 617722) = . 

We are therefore left with different possibilities, that we study in the next sec- 
tion. Note that, although the denominators in eq. flT^ vanish when m\ = m\, 
we will see in the next section that no singularity arises, since the numerator 
then also vanishes, because ci = 62. 

3 Solutions of the Gap Equations - Dynami- 
cal Fermion Masses and Mixing 

We now detail the different solutions to the gap equations flT5|) . The 
trivial solution corresponds to the situation where no dynamical mass is 
generated, tt^i = m2 = /^ = 0, which is of no interest to us here. In what 
follows we focus on situations, in which fermion masses are generated with 
the constraints (pPj) satisfied. 

3.1 The case mi = m2 = and /i 7^ 

In this case, the eigen masses are 

A± = ±/x , (21) 
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and the mass eigenstates are 

V'± = -^(V^2±^i) , (22) 

such that the mixing angle (TTOl) is 6' = — 7r/4, in our conventions. 
This case does not include a mass hierarchy, hence there are no oscillations 
(TTT]) among the fermion flavours either, since the energy eigenvalues E = 
\/'p^ + /i^ are the same. 

Among the four gap equations (fT5!) in this case, only the third is not 
trivial, and leads to 

(eie2-e2)(Ji-^2j^) . (23) 



4 + C 
Since A^ = /i^, the expressions (fl8|) lead to 



and we obtain 



h 


r^ 


-1 
167r2 


/ 

■ 1 + In 1 

V 


^M2 




h 


— 


-1 


1 






167r2 


/i^' 




In 


,M2 




-167r2 






(4 + C)(e: 


[62 - 


e^) 



(24) 
(25) 

(26) 



We note that this expression has a meaning only if 6162 > e^, otherwise 
/x^ > M^. Assuming this constraint on the couplings, we finally obtain 

/i-Mexp( ~^''' -) . (27) 

V(4 + C)(eie2-e2)y 

3.2 The case e2mi + eim2 = and m^ 7^ m| 



In this situation, the first equation f lT9|) leads to Ji = 0. The expression 
( IT8|) for Ji leads then to 

A^ = A^ = exp(-l)M2 , (28) 

which is not physical, because the dynamical masses are then necessarily of 
the order M. This possibility is therefore disregarded, since we will eventually 
take the limit M — )■ 00. 
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3.3 The case rrii = — m2 7^ 

In order to have mi = —m2 = m, it can be seen from eqs.f llSp that 
necessarily ei = 62, such that both constraints fl20|) are satisfied. Also, 
eqs. flT5|l are equivalent to 

(e'-e')[h-{fi' + m')h], (29) 



4 + C 
and A^ = m? + /i^, such that we find 

m^ + ^^ = M^exp( ^^^-y_^,^ ), (30) 

which has a meaning only if e^ > e^. This condition allows one to take 
the limit e — )> without affecting the mass eigenvalues or mixing angles 
(see below). This is important, because, as already mentioned, a non-zero 
flavour-mixing coupling e might lead to dynamical generation of vector boson 
masses [T9l [20] , thereby spoiling their nature as regulator flelds. 

We stress here that we cannot determine m and /^ independently, and the 
eigen masses are 



A± = ±^m? + /i2 . (31) 

The mass eigenstates are 

^± = ^ f^i + ^rT^==.A ' (32) 

m ± y m^ + /x / 




where 

2 _ 2m2 + 2/^2 ± 2mv/m2 + /i2 



Nl = ' "" -"-V- ' ^ ^ (33) 

2771^ + /i^ ± ^ra^m?' + /i^ 

and the mixing angle ^ (fTOj) is given by 

tan^ = 7^^ . (34) 

In order to flx the mixing angle, one would need an additional ingredient, 
since the present model does not flx /i, but only m} + [j?. 

Again, there is no mass hierarchy due to ( ISTl) in this case, the energy 
eigenvalues are the same, and so no oscillations ( ITT]) among fermion flavours. 
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3.4 The case rrii = 1712 ^ 0: Dynamical Flavour Oscilla- 
tions 

We find here from eqs. flT^ that necessarily ei = 62, e = and /x^ = m^. 
we have then 

H^ = 71111712 = m^ and h = ,. ,. „ , (35) 

(4 + Qe^ 

where e = ei = 62- We have A^ = and A^ = 4m^, such that the expressions 
(IHD and dSSD for h lead to 

' ^n(^^)^,^±^, (36) 



167r2 \M'^ J (4 + C)e2 
and the common dynamical mass is finally 

which, as expected, is not perturbative in e. In this situation, the mass 
matrix has identical elements, and has the eigenvalues 

X^ = 2m = Mexp(- ^^_^ ] , A_ = , (38) 

and the corresponding mass eigenstates are also given by eq. fl2^ . The mixing 
angle ( ITO!) is ^ = =F 7r/4, depending on the sign of /i = ± m, respectively. 

In this case, one of the fermions is massless, and the other massive, with 
mass 2771. There is a mass hierarchy and thus oscillations (TTTT) among the 
fermion fiavours in this case. We note that because of the constraints (l2Up . 
this is the only case in the model ([3]) where non trivial oscillations among 
fermion fiavours take place. As we have seen above, in this case necessarily 
the fiavour-mixing gauge couplings e — )■ 0, so one does not have to worry 
about dynamical generation of gauge boson masses, and thus the latter play 
a consistent role as regulator fields. 

3.5 The case e = 0, /i = 

This is a straightforward generalisation of the original model of |1], which 
involved one fermion, to the two fermion-fiavour case with no mixing at all. 
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This case can be divided into two situations: i) rrii ^ and m2 7^ 0, and ii) 
nil = or 7712 = 0. As we shall discuss in section HJ these may be relevant for 
Majorana neutrinos in the standard model or extensions thereof, involving 
right-handed neutrinos, respectively. 

3.5.1 i) mi 7^ and m2 7^ 

In this first situation, we obtain from (|9]) that the two eigenvalues of the 
mass matrix are 

7ni = Mexp( ~^'[ j ,^ = 1,2, (39) 

so the dynamically generated mass matrix is diagonal with masses TUi among 
the two flavours. Hence there is no mixing ( ITOl) or oscillations (11 II) between 
the flavours ipi, i = 1,2, in this case. 

3.5.2 ii) m-i = or m2 = 

In this case, we observe from the system of equations flTSj) that there is 
also a consistent solution, with either ttii = with 7712 7^ or 7n2 = and 
7711 7^ 0. The two cases are symmetric. For reasons that will become clear 
from our discussion on neutrinos in section HI we may concentrate for brevity 
in the former case, i.e. 7ni = 0. In that case, the solution of equations (fT5!) 
yields 



while from (IT8|) and the deflnitions fITTj) we obtain that in this case A^ = 0, 
while A"^ = 7712 7^ 0, and thus from (jlQj) we have 

|^,| = M==Afexp(-^j|^) . (41) 

Note that, although I2 diverges logarithmically as /i — > 0, it enters the gap 
equations (fT5|) only in the combination {fi'^ — 7ni7n2) I2 = IJ^Iii which vanishes 
in this limit. 

The mass eigenvalues are in this case, A_ = 0, and A+ = m2 7^ given 
by (141 p . The mixing angle 9 is though vanishing and thus there are no 
oscillations between the states. 
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3.6 Lorentz symmetric limit 

In order to recover Lorentz invariance, we finally take the simultaneous 
limits 

M -)■ oo and ei, 62, e -^ , (42) 

in such a way that the dynamical masses are finite, and we denote the corre- 
sponding "renormalized" mass matrix by Mr. This procedure is independent 
of the gauge parameter (, and the resulting fermion mass is set to any desired 
value. In this limit, the gauge field decouples from fermions, and the only 
finite effect from Lorentz violation in the original model is the presence of 
finite dynamical masses for fermions. 

We now check this statement by demonstrating that the fermion disper- 
sion relations are relativistic in the limit fH21). We focus here for concreteness 



on the solution described in subsection 13. 4[ with fx = +m, but clearly the 
same conclusion holds for all the other solutions given in section |3l Because 
one of the eigen masses vanishes, which leads to one-loop infrared (IR) di- 
vergence, we consider mi = m2 = m and m — fi = m5 with 5 << L As will 
be seen, after the limit (142 p is taken, the fermion self energy won't depend 
on (5, such that the limit 5 — )■ will not introduce any IR divergence. We 
calculate in Appendix B the one-loop fermion self energy, where we use the 
Feynman gauge since the limit (H2l) is gauge independent. To lowest order 
in momentum, we find then 

^={^h fh)u^l'-(fr fl^Ap-l-M, (43) 

\^off ^diagj K^off ^diag/ 

where {u},p) is the external 4-momentum and 



e^/1 1.^ 1. . ,/m 



Z" = ^ In 2 + - In 5 + In— (44) 

'^"'^ 87r2 V4 2 2 \mJJ ^ ' 

•^'"3 87r2 V 12 2 2 \m) J 



e" 



Z'off = Zl,f = ^ln2-\n6). 
As expected, because of Lorentz-symmetry violation, Z^^^^ 7^ ^diag' but since 



^'nMj = -2-^ («) 
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the limit (HZD leads to 



E ^ --(u;7°-p-7)l-Mfl . (46) 

Therefore the dispersion relations are relativistic, since time and space deriva- 
tives are dressed with the same corrections in the limit ( l42l) . These correc- 
tions can be absorbed in a fermion field redefinition, so that we are left with 
two free relativistic fermion fiavours oscillating. 

3.7 Energetics Arguments 

Among the different possibilities to generate masses dynamically, one can 
question the preference for the system to have non-vanishing masses, rather 
than no dynamical mass generated. We give here an energetics argument 
supporting the choice of non- vanishing dynamical masses [20]. This argument 
is based on the Feynman-Hellmann theorem [22], which states that, if there 
is a ground state |\1/a) of a system with Hamiltonian that depends on a 
parameter A, then for the energy E of this ground state we have: 

^=(^.1^1^.), (47) 

where H is the Hamiltonian operator of the system. In our situation, the 
parameter A can be chosen to be A = M~^, such that 



^ m(0| y^^x^ {F,,^F^-)j, |0)m , A = M-^ (48) 



dE 
'dX'' '4 

where the index E denotes Euclidean formalism, as a result of the fact that 
the Hamiltonian of the system is identified with minus the effective Euclidean 
action. One should expect that the Lorentz-violating nature of the vacuum 
|0)m implies in general the non vanishing of the right-hand-side, implying a 
dependence of the vacuum energy on the dynamically generated mass. Using 
the cyclic Bianchi identity for the gauge bosons field strengths, 

5[M^.p] = 0, (49) 

with the symbol [. . . ] denoting symmetrisation of the appropriate indices, we 
obtain 



dE _ 1 
9A " "4 ^'^ 



cf'x^ {F^,d,[d^F'' + d'^F'^])^ |0)m . (50) 
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Integrating by part and assuming that the fields decay away at space-time 
infinity, one may write eq. fl48p in the form: 

^ = +1 MiOlJd'x"^ {d^F,^d,F-)^ |0)m (51) 

We write then the equations of motion for the vector fields, from the La- 
grangian ([3]) where we neglect the operator A/M^, and we obtain: 

^ = +1 m(0| y"rf^x^ [{jy + J- J- JkdoF'')^ |0)m , (52) 

where the current is J^ = \[f7^r\E' and the Euclidean formalism used in fl52|) . 
In the framework of the LIV model studied here, one might face a situation 
where non-trivial condensates of the covariant square of the stationary four- 
current J^ are observed in the (rotationally invariant) vacuum. For such 
stationary currents, where d^F^^ = 0, we have then from eq. (!52|) : 

-q^ = ^m{0\J d'x^{J^J,)^\0),,>0. (53) 

This implies that the vacuum energy E in this case is a monotonically de- 
creasing function of M^, so that the energy goes to its minimum in the 
Lorentz symmetric limit M — ;• oo we are interested in. 

The above argument in favour of the stability of the Lorentz Invariant 
Limit ([2]) can be turned into an argument in favour also of the dynamical 
fermion-mass generation as follows: in a finite M < oo situation, the gauge 
coupling e is considered as an independent quantity from M, and thus, in 
view of dl]), the LIV mass scale is proportional to the fermion mass m > 
(absolute value if ?72 < 0). In this sense, one obtains from fl53j) . 

dE dM dX dE 1 



M 



{0\Jd^x''{J>'J,)^\0)M<0. (54) 



dm dm dM dX mM"^ 

Thus, the energy of the vacuum for any finite value of M is also a mono- 
tonically decreasing function of the fermion mass. In the Lorentz-symmetric 
limit dl]), the energy exhibits a plateaux, as far as its dependence on the 
finite m > is concerned, i.e. dE/dm = 0, but its value is lower than the 
case where ?n = 0. 

We stress, however, that the above arguments rely on the formation of 
condensates for the covariant square of the current. The latter property is 
at present a conjecture, and its proof goes far beyond our considerations in 
this article. 
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4 Extension to Chiral Majorana Neutrinos 

Above we considered Dirac non-chiral ferinions. However, if we wish to 
present the above-described dynamical mass generation scenario as a viable 
alternatives to standard seesaw mechanisms for neutrinos, and explain the 
neutrino oscillations as a dynamical phenomenon, then we should extend the 
above considerations to the case where the fermions are chiral and Majorana 
(as most likely is the case realised in nature). 

Below we shall consider two separate cases. The first is the one in which 
the fermions correspond to Majorana mass eigenstates obtained from the 
left-handed flavour neutrino physical fields of the standard model, while the 
second case involves sterile right-handed neutrinos as in seesaw extensions of 
the standard model. 

4.1 Left-handed Neutrino Majorana Mass Generation 

For instructive purposes it is useful first to review some basic formal- 
ism. According to the standard theory [23] a (Majorana (M)) mass term for 
neutrinos, which involves only left-handed fields, reads 

C'^ = -IvzM'' {ulT + \^.c. . (55) 

where the normalisation of 1/2 will be understood in what follows. In the 
one generation case we focus upon here M^ is a c-number (In case of many 
generations, z/^, i = e, fi,T then M^ is a symmetric 3x3 matrix, as can be 
seen easily). The (physical) Majorana field u^^ , involving both chiralities, is 
defined as 

u'' = UL + {vlY (56) 

and is always an eigenstate of the mass, that is when expressed in terms of 
it the Mass matrix is diagonal 

£M ^ _ l^M^ ^M ^J_J2m,17,u,, (57) 

with TTii the mass eigenvalues. It satisfies the Majorana condition 

(u^^^Y = u^ , (58) 
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which imphes that a Majorana field is its own antiparticle. 

The kinetic (Dirac) term of the Lagrangian with respect the left-handed 
ul fields, when expressed in terms of the Majorana mass eigenstate fields Vi 
reads (up to an irrelevant total derivative): 



^Ckin = Vhi^i^L + h.c. = ^ T^^jif^^j + h.c. (59) 

j 
with Uj = {ujY = Uj L + {i^j lY , 



with the suffix j denoting mass eigenstate fields. In view of the extra | nor- 
malisation of the kinetic terms of the Majorana fields then, it is customary to 
define the corresponding mass terms with the same normalisation, as we have 
done above, compared to other Dirac fields one encounters in the standard 
model. 

In our approach we consider the coupling of a doublet of (mass eigenstate) 
Majorana fields to the regulator U(l) gauge field A^ in the case discussed 
in subsection 13.5.11 The fact that a Majorana field contains both chiralities 
allows for a straightforward extension of the Dirac case discussed in previous 
sections to the current situation. In this way, we are able to generate dynam- 
ically different mass eigenvalues ( 139|) for the two species, without mixing, as 
implied by the solution of that section. This is a consistent way of discussing 
the dynamical appearance of a Majorana mass for left-handed neutrinos of 
the standard model. Non-trivial mixing of flavour neutrinos, coupled to the 
physical SU{2)l gauge flelds of the standard model, can then be obtained 
in the case where the mass eigenvalues are different. In order to recover 
the Lorentz symmetric limit in this case, we need to take simultaneously 
Ci, 62 — )■ in such a way that their ratio is flxed to the phenomenologically 
desired value. It is important that in this approach we started from Majo- 
rana mass eigenstates coupled to the regulator gauge flelds, with no mixing. 
The latter is obtained when one expresses the Majorana mass eigenstates in 
terms of the flavour neutrino eigenstates, which appear in nature. 

4.2 Extensions of the standard model with right-handed 
(sterile) neutrinos 

When there are right-handed (sterile) neutrino components present, vr, 
one can deflne two kinds of mass terms, Majorana ( 155|) (M) and Dirac (D). 
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The most general mass term, then, reads 

/:''+^ = -]^L Mf {vlY - Vl M^ un - \vr M^ {urY + h.c. , (60) 

where the mass matrices Mffj^ and M^ are in general different. 

We consider below the mixed mass terms (1601) in the one generation case, 
of relevance to our models discussed in this work. In this case we may as- 
sembly the left-handed neutrino fields and the conjugate of the right-handed 
one into a left-handed doublet field 

in which case the mass term (I60p can be written in terms of a 2 x 2 mass 
matrix (6 x 6 in the case of three generations): 

where for the sake of generality we expressed here the mass matrix as a 
matrix with flavour components as well. For a single generation of neutrinos, 
we consider below, the elements of the above (2 x 2 in this case) matrix 
are c-numbers. For our toy purposes here we assume no CP violation in the 
lepton sector [23] . 

The matrix M^"*"^ can be diagonalised by a Hermitean matrix U: 

with 1231: 



U = Or^'l\ 0=(''''\'^^^ ^ (63) 

y— smt' cos a J 

where the matrix rj has eigenvalues rji = ±1, which are related to the so- 
called CP parity of the Majorana neutrinos [23], and stem from the fact that 
the mass eigenvalues can be positive or negative 

<,2 = \ {Mr + Ml) t \ ^J{MR-MLf + ^Ml (64) 

so one can rewrite them as m.j = |m.j|?7j = frii ?7j, rji = ±1, i = 1,2. The mixing 
angle 9 being such that: 

-— , ^"-^^^ ta„2»=^M^. (66) 

^{Ma - Ml,)'' + 4Ml Mh-Ml 



The Majorana fields, involving both chiralities, are then defined in terms of 

U as 

ly^ = U^nL + (U^ ulY = h) , ly'r = ly^, ^ = 1, 2 . (66) 

These are the mass eigenstate fields with masses mi^2 dHlD- 

The original chiral (left-handed) neutrinos, appearing in the Lagrangian 
(155|) are related therefore to these mass eigenstates as follows: 

ul = cos 6y/rj^ui L + sin 6y/rj^U2 l 
{i^rT = -sm9y/rjluiL + cos9y/rj^U2L (67) 

In the standard seesaw scenarios [2], there are no masses for the left handed 
fields. Ml = 0, and the right-handed neutrino (sterile) Majorana masses are 
assumed much heavier than the Dirac masses , M/j ^ Md, the latter being 
given by means of a Higgs mechanism by, e.g., Yukawa coupling terms of the 
form 

(I)^Vli^r + h.c., (68) 



where cj)'^ = ia2(p* is the dual of the Higgs doublet. In this limit, from ([& 
(!65|) the mass eigenstates generated are of the form mi ^ jf- <^ Mn and 

^^^ 1 an 

Mr 



m2 ^ Mji :$> Mf), while the mixing angle 9 ~ ^ ^ 1, and also rji = —1 



7^2 = 1, hence from f l67|l we do obtain: 



, Mn 

{vrT - -iYJ^J^iL + Jy2L ■ (69) 

The purpose of the remainder of this section is to adopt the previous pro- 
cedure and generate dynamically masses for the Majorana fields by coupling 
them to gauge fields. To this end we view one of the fiavours as a right-handed 

sterile neutrino, Nji = ^ ( I-I-75 ) A^ and the other fiavour ^l = |( 1— 75 ) i^, as 
an active neutrino of the standard model. Here, N, ijj are non-chiral spinors, 
which in the case of neutrino may be taken to be Majorana. This would be 
a toy (with one active and one sterile neutrino) version of the minimal (non 
supersymmetric) extension of the standard model of ref. [23], termed z/MSM. 
In this case we do not avoid right-handed neutrinos but we use the dynamical 



19 



mass generation mechanism presented here to give masses to them. In this 
case the Lagrangian (^ is replaced by: 

^ = -\F,ua-^)F"' + N{^<^-e^4)^-[l + ^,)N (70) 

+ V^(z^ - 624) - (l - 75) ^ , 

Notice that in this case, due to the opposite chirahties of the two spinor fields, 
the off diagonal fiavour mixing gauge couplings e are irrelevant because the 
corresponding terms vanish identically. 

According to our general discussion on combined Dirac and Majorana 
masses above, we may express this Lagrangian in terms of Majorana fields 
and view the initially massless ip and A^ as the Majorana field doublet u'^ 
fl66|) which then couples to the vector fields. Dynamically generated mixing 
of the two should involve a small mixing angle in phenomenologically realistic 
situations in view of the discussion above (jHSD- 

Unfortunately, in our single gauge field toy models considered here, the 
only solution from the cases discussed in section [3] that can be carried over 
to the case of Majorana neutrinos is the one discussed in subsection 13.5.21 
In this case, the quantity mi = can be identified with the left-handed 
Majorana mass Ml = 0, which in the usual seesaw models is assumed zero. 
The dynamically generated mass eigenvalue A+ = m2 given by eq. fj^T]) is then 
identified with the heavy right-handed Majorana mass, Mr, but there is no 
non-trivial Dirac mass /i, since the latter vanishes in the dynamical solution. 
The latter can be generated through the usual Yukawa couphngs (jSHD with 
the Higgs field, which upon acquiring a vacuum expectation value via the 
Higgs mechanism would generate a Dirac mass term. 

5 Conclusions and Outlook 

In this work we have considered the coupling of fiavoured fermion fields 
to LIV vector gauge bosons, with Lorentz Invariance being violated in the 
gauge sector at a mass scale M and studied the limiting case where the gauge 
couplings go to zero, while the LIV Mass scale M — )■ 00 simultaneously, in 
such a way that the Schwinger-Dyson dynamically generated fermion masses 
remain finite. No vector boson mass is generated due to an appropriate ar- 
rangement of the couplings. In this way, the LIV vector bosons are viewed 
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as regulator fields, with the only remnant of the LIV the dynamical fermion 
mass. Unfortunately, the dynamical equations are sufficiently restricted so 
as to allow only one case where oscillation among fermion flavours is al- 
lowed and in this case one of the fermion mass eigenstates is massless, while 
the other is massive. The mixing angle is necessarily maximal in this case 
9 = ±7r/4. One may hope that extension to a third flavour may lead to 
more phenomenologically realistic situations with arbitrary mixing and mass 
generated for all flavours. 

Another possibility towards this result might be the inclusion of more than 
one regulator vector flelds, along the lines of [20], where, however, not only 
a mass hierarchy is generated between the fermions, with non trivial masses, 
but also one of the gauge bosons acquires a mass. In our case of regulators, 
unfortunately, this last mass would also be kept flnite, but probably this 
would not be a problem, since the massive vector fleld decouples from the 
Lagrangian of the fermions in the zero gauge coupling "relativistic limit" d?]). 
We hope to come back to this case in a future work. 

Another aspect of our work, which was also the original motivation, is the 
one in which this method applies to chiral neutrinos of the standard model, in 
an attempt to discuss neutrino mass generation independently of the seesaw 
mechanism. We have discussed two scenarios in this respect. 

In the flrst, we avoided the inclusion of sterile neutrinos altogether. In this 
case the two flavours considered above have been viewed as corresponding to 
Majorana mass eigenstates of two left-handed neutrino flavours, interacting 
with a LIV regulator gauge fleld with vanishing couplings. It was demon- 
strated that different mass eigenstates could be obtained in the Lorentz sym- 
metric limit, which then leads to standard oscillations among the physical 
neutrino flavours coupled to the SU{2)l gauge flelds of the toy standard 
model involving only two flavours. Extension to the physical case of three 
generations, including CP violation in the lepton sector, will constitute the 
subject of a forthcoming publication. 

The second scenario, involved an extension of the model (j3]) to a toy 
version of the i/MSM model of [21], in which one of fermion flavours of (E]), 
say ipi, represented a right-handed neutrino fleld, and the other flavour 1JJ2 
a left-handed active neutrino of the standard model. In such a case our 
aim was to generate dynamically a mass hierarchy between active and right- 
handed (possibly sterile) neutrinos of the type needed in phenomenological 
approaches to dark matter, where a keV sterile neutrino may play the role of 
a dark matter fleld, consistently with current astrophysical and cosmological 
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data [23] . In the context of our framework, we can only generate dynamically 
the right-handed neutrino mass, but not a Dirac mass term. It is interesting 
that the absence of a left-handed Majorana mass (standard assumption in 
seesaw models) appears naturally in our models. A Dirac mass then, coupling 
left (active) and right-handed(sterile) components can be generated by the 
standard Higgs mechanism of the standard model. 
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A Appendix: gap equations 

The aim of this appendix is to present the main steps to obtain (fT5l) from 
the Schwinger-Dyson equation (fT4l) which is rewritten below: 

G-'-S~'= j D^^T^GrY . (71) 

Jp 



The first step that we take is to commute the first r and 7'^ in fjTTI) . so 
that in the middle of the integrand we have a matrix product given by 

tGt = X h -'"] f^-"^' ^ "l ('} -''] (72) 

^ fel{f - nii) + e^{f - mi) -Y + Ai(eie2 - e^) 

F + /i(eie2 - e^) e^ {f - m2) + el{f - mi\ 
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where 



_ . p^ + f{mi + 1712) + mim2 - fx"^ , . 

(p2 — rnl)(p'^ — ml) — 2^'^{p^ + 71111712) + //^ 

p"^ + p{mi + 7722) + 77117712 — /i^ 

= z— 



y = ie[ei{f - m2) + e2{f - mi)] , 

with A^ defined as in ( fTTl) . If we identify individually each matrix element 
in the Schwinger-Dyson equation (17T1) . we obtain for the Mn element 

zmi = I D^,^>'X[el{p-m2) + e\f-mi)]Y (74) 

(4 + C) p^(efmi + e^m2) + (/i^ — mim2)(efm2 + e^mi) 



The last equation can be written as 

'^l /.2_ I ^2_ \T I /,,2 _ _ ^/'^2_ I ^2„ 



(e^mi + e m2)/i + (/i - mim2){eim2 + e mi)l2, (75) 



4 + C 
where 

1 1 

/p 1 +p2/M2(j92-A^)(p2 _ A^) 



h = -I I -v-^nTT^— — T2\n:2 — jr^ (^6) 



/2 = -I 



,/^l+p'/M^)p\p^-Al){j)^-A\) ■ 
The Wick rotation p^ -^ iu leads to 

J ^ 1 r fdp r- du 

' 47r3yo l+p2/M2y_^(c^2+p2 + ^2^)(^2+p2 + ^2_) I ) 

1 p p^rfp p du 



The integrand of Ji which depends on u only can be written 



1 



^2-^2 



1 



1 1 



a;2+p2 (^2 + p2 + ^^ 
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(78) 



and similarly, the integrand of I2 which depends on u only can be expressed 

as 

1 1 



Al - Al [Al \uj^ + f {u^+f + Al 

1/1 1 



Al \lo^ +p^ (a;2 + ^ + A^ 



(79) 



Finally, substituting (!78|) and (179!) into (1771) leads to the first equation of 

(USD. 

Furthermore, due to the symmetry of our model, the second equation of IHM 
is obtained from the first one by exchanging mi and m2. 
Finally, the left-hand side of ( 17T|) is symmetric with non-diagonal elements 
given by i/x, therefore, the non-diagonal elements of the right-hand side must 
also be equal. However, looking at ( 172|) . we realize that it is only possible if 
the terms related with Y vanish. So, the non-diagonal elements give us the 
following equations 

/■ (4 + C) . 2N P^ + mim2-yu2 

'^ = i (1 + f/M^^'"'' ~ ' \\P- - AlW - Al) ' ^'°^ 



= ID^.^^'XYY 

(4 + C) p^(eimi -I- 627712) + (/i^ - mim2){eim2 + e2'mi] 



[l+f/M^) p^{p^ - A'i){p^ - Al) 

where using eq. (176l) . leads to the last two equations (1T5|) . 

B Appendix: one-loop fermion self energy 

We calculate here, in the Feynman gauge, the fermion wave function 
renormalization for the case {ei = 62 and e = 0}. In order to avoid IR 
divergences obtained in the one-loop calculation for mi = m2 = /i, because 
one of the eigen masses vanishes, we consider here the situation mi = m.2 = 
m ^ fi. The fermion propagator is then given by 

G(p) = ^ P' + ^^P + ^'-f^' (p-m /i \ 
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We obtain the fermion wave function renormalization by differentiating the 
fermion self-energy with respect to the external momentum and then, set 
it to zero. Since the fermion propagator f lHT]) has two independent flavour 
components, we consider the one-loop diagonal self energy S)^^^ and the 
one-loop off-diagonal part SJ^A, where 

^Z(^^P) (82) 



-le 



^ '' d^k f l^l^J.[{p — kY — {w? — iJ?)]m 



(27r)4 j 1 + piM-i I k^[{p - kf - (m + /i)2][(p - kf - (m - /i)2] 
Y{f - ^HUp - kf - K + /i2)] 



+ 



k"^ [{p — ky — {m + jiY] [{p — kY — {m — jiY] 

off 



^offi^^P) 



-le^ f d^k 7'^7/,[(p-A:)2 + m2-/i2]/i + 2m/i7'^(^-^)7^ 



{2-kY J 1 + k'^/M'^ k'^ [{p -ky -{m + iiy] [{p - kf - {m - jif] 

Differentiating these terms with respect to Pp and then setting w = and 
p = 0, we find 



d^^diag\ _ ie^ f d^k f A;V - (m2 + ^2)^p ^ 2A;^^ 



dpp 'P=o Svr^ y I j^ p I M"^ \k'^[k'^ - {m + ^YWk"^ - {m - ^Y] 

_ AkP^k^ - 8kPi(:k\m^ + fi^) + 4kf^{m^ + fi^ ' 
k^[k^ - (m + /i)2]2[A;2 - (m - ;u)2]2 

^^off I ifime^ f d'^k f —7'' 



(83) 



dpp 'P=o Att^ J i + ^2/7Vf2 [A;2[A;2- (m + /i)2][A;2- (m-^)2] 

AkP^k"^ - AkPl^{m'^ + /i2) 



+ 



fc2[/c2 _ (^^ _(_ ^)2j2[/.2 _ (-^ _ ^pj5 



(84) 



We write then 



^diag = -^ + Zdiag^l - Z^^^^p ■ -f 

Sl;^/ = -/i + ^V7°-^oV-7, (85) 

and since we are interested in the limit /i — )■ m, we write m — fi = mS, with 
S « 1 and approximate m+jj ^ 2m. In terms of new variables x = yk'^/m, 
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y = ko/m, A = m/M <^ 1 and after a Wick rotation, we obtain 



7° 

diag 



27r3 



x'^dx 

1 + A2x2 



rfy 



-(x2 + 2/2) - 2y2 - 2 



_(x2 + 1/2) (2;2 + y2 _|_ 4) (^^2 + y2 _,_ ^2) 



2/'^2 



4y^{x 



y 



2\2 



16r(a;^ + ?/^) + 16r 



(X2 + y2)(a;2 + ^2 + 4)2(3.2 + ^2 + ^2)2 



(86) 



7° 



and 

diag 



-2 /"^ a;2dx 



Tpo 







1 + A2x2 



dy 



_(X2 + 1/2) (a;2 + ^2 + 4) (3.2 + ^2 + 52) 



%2(x2+?/2) _^3^2 



(a;2 + 2/2)(^2;2 + y2 _^ 4)2(a;2 + y^ + 5'^y _ 
"°° x^dx 



27r3 







1 + A2x2 



dy 



-(x2 + 2/2) 



2x2/3 



4 a;2(a;2 + 2/2)2 .^43,2(^2 _j_ ^2) _^ 4^2 



7I 



7J-U 



3 (x2 + 2/2) (x2 + 2/2 + 4)2(x2 + y^ + 52)2 
dy 



1 + A2x2 



_(X2 + 2/2) (x2 + 2/2 + 4)(x2 + 2/2 + 52) 

(87) 

1 



X2 + 2/2) (x2 + 2/2 + 4)(x2 + 2/2 + ^2) 



4 



x^(x^ + 2/ ) + 2x^ 



3 (X2 + 2/2)(x2 + 2/2 + 4)2(x2 + 2/2 + ^2)2 

Finally, after solving these integrals, we find 



7° 

diag 






^vr2Vi-^1^2 + lln5 + lnA 



(In 2 -In (5) 



167r2 



and 



7I 

diag 



7I 



87r2 

e2 

167r^ 



1 1 1 

In2 + -In^ + lnA 

12 2 2 



(In 2 -In 5) 



(89) 
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